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A simple molecular model of chiral molecules is presented in this paper : the chiral dipolar 
hard sphere model. The discriminatory interaction between enantiomers is represented by 
electrostatic (or magnetic) dipolos-dipoles interactions : short ranged steric repulsion are 
represented by hard sphere potential and, in each molecule, two point dipoles are located 
inside the sphere. The model is described in detail and some of its elementary properties 
are given ; in particular, it is shown that the that the knowledge of only three multipole 
spherical components (namely : Qio, Q21 and Q22) allows to compute all multipole spherical 
components of the model. Despite, the simplicity of the model, it is shown also that the energy 
landscape of the interaction between two enantiomers is quite rich, this renders systems of 
chiral dipolar hard sphere very interesting and complicated to study. Few preliminary Monte 
Carlo simulation results are also reported in the paper. Last, but not least, this paper is 
dedicated to Jean-Jacques Weis. 

1. Introduction 

Chirality is used to describe a lack of symmetry in compounds : it refers to ob- 
jects, molecules or structures that are not superposable on their mirror image. The 
separation of the two optical active forms of the tartaric acid, based on the obser- 
vation of the symmetry of crystals by Pasteur in 1847, has been a major step in 
the understanding of optical properties of crystals and liquids [ij]. Later, in 1884, 
Kelvin introduced the term chirality, but it has been extensively and systemati- 
cally used only eighty years later with the introduction of the Cahn-Ingold-Prelog 
priority rules [4]. Since Arago, Biot and Pasteur, the ideas and concepts on optical 
activity of substances and on chirality have had and still have a lot of applications 
in chemistry, biology and physics. 

In chemistry and biochemistry, chirality is an important tool in stereochemistry 
and for the determination of the functions of molecules. The interactions between 
biological molecules are frequently a consequence of the chiral interactions between 
them. For instance, as first noted by Pasteur, the two forms of asparagine {dextro 
and levo) tasted differently ; it was found later that this is a frequent property 
the two optical active forms of amino acids, because of the chiral structure of the 
taste bud receptor site (several others interesting examples may be found in [1]). 
The helical structures of DNA, RNA, proteins, etc. are also important to permit 
biological recognition and interactions between biological structures [1, [H| ■ 
In condensed matter physics many properties depend on the chirality of consti- 
tuants or structures. More specifically, in liquid crystal physics, many phases as 
cholesteric, some smectics phases {S\, Sq, etc.) and blue phases possess an optical 
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rotatory power [1, Q, @, S ; in most cases, a helical arrangement of the 

molecules is responsible for this property. Obviously, these properties are of a par- 
ticular interest for technological applications. It is not necessary for the molecules 
to be chiral to observe chiral liquid crystal phases 0, [13, [il ; however, some par- 
ticular phases, such as the ferroelectric smectic-C* phase, may be observed only 
when some chiral molecules are present in the sample (see section 5.10 of ref.[9[] 
and refs.[il, 

To describe the optical activity of liquids and solids, and the discriminatory inter- 
action between chiral molecules several molecular models have been built. Among 
the most known are Kuhn models that are based on coupled-oscillators arranged on 
an asymmetric tetrahedron [l^. the optical activity tensor or gyration polar- 
izability may be computed in these models 13, U, 22]. Models based on chirality 



of helix have also been used [21|, |23j, |2j] . These helical models are important for 
the study of the pro perties of carbon nanotubes j25l . |26| . 1271 ] and biological macro- 
molecules 0, m, [2^ . 

For the study of chirality in fluids and solids, the discriminatory interactions be- 
tween enantiomers are important too 0]. The difference in the interaction between 
enantiomeric molecules may have several origins, it may stems from steric, electro- 
static, magnetic or dispersion interactions ja. l3ll. |32|]. These chirodiastaltic interac- 
tions [30(] are responsible for some differences in thermodynamical and structural 



properties between enantiopure systems (homochiral systems - systems that are 
composed with exclusively one enantiomer) and racemic systems (systems com- 
posed in equal proportions of the two enantiomers). For instance, the melting 
point of crystals of the enantiomeric pure (-l-)-tartaric acid is 170°C, while, for 
the racemic crystals, the melting temperature is 204 — 6°C. Some differences are 
also reported for boiling points, but since the chirodiastaltic interactions are quite 
small in general, these difference in boiling points are small and do not occur sys- 
tematically. In liquid crystals, the phase diagrams depend also on the fractional 
concentrations of enantiomers [s^ ]. 

Several simple molecular models have been proposed and used to study the ther- 
modynamical and structural properties of fluids and solids of chiral molecules in 
computer simulations and integral equation theories. In these simple molecular 
models, the chirodiastaltic interactions are represented either via steric interac- 
tions [H, 35, m, 37, 38] or electrostatic interactions [1^, [13, 41, 42] ; a chiral 
Gay-Berne model has also been built [i^]. In the present work, we describe an- 
other simple model of chiral molecules, the chiral dipolar hard sphere model. In 
this model, the chirodiastaltic interactions is represented by electrostatic (or mag- 
netic) dipoles-dipoles interactions, it consists of hard sphere with two point dipoles 
located inside the molecule. As shown by Craig and coworkers, electrostatic interac- 
tions between enantiomers may be discriminatory at most if the dipole-quadrupole 
interaction between two enantiomers differs [3, 0, [4^ ; this implies that if the 
quadrupole moment of the two form of enantiomers differs then the interaction 
may be discriminatory. Therefore, a chiral model based on electrostatic interac- 
tions must have at least four point electric charges arranged on an asymmetric 
tetrahedron geometry, as the models in refs. non-colinear point 

dipoles 461] or one quadrupole with no planar symmetry 44]. 
The present paper is organised as follow. In section 2, we describe the model of 
chiral dipolar hard sphere in detail and give some of its elementary properties. 
More precisely, in section 2.1, we define the model, the nomenclature, its pertinent 
parameters and we compute all its multipole moments ; it is shown that the knowl- 
edge of only three multipole spherical components (namely : Qio, Q21 and Q22) 
allows to compute all multipole spherical components of the model. The knowledge 
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of all multipole moments is important in the definition of the model and in the de- 
termination of its electrical properties and also, in view of the determination of the 
optical activity of the model, for the computation of multipole-multipole polariz- 
abilities [s, 19, 44, 45, 47, 48, 49, 5^. In section 2.2, we compute interaction energies 
between two enantiomers in some particular configurations ; it is shown that, de- 
spite the simplicity of the molecular model, the energy landscape of the interaction 
between two enantiomers is quite complicated. In section 3, we report some prelim- 
inary results of Monte Carlo simulations performed in the canonical (NVT) and 
isobaric (NPT) ensembles ; in these computations, all dipole-dipole interactions 
between molecules are explicitly taken into account by using the Ewald method 
for dipolar interactions [51[. The paper ends with a discussion on the perspectives 
that this model offers to the study of chiral fluids and solids. For completeness, in 
appendix A, we give general formulas for the computation of multipole moments 
for any dipole distribution. 

Last, but not least, this paper is dedicaded to my colleague Jean-Jacques Weis. 



2. Elementary properties of chiral dipolar hard sphere model. 

In this section, we describe the chiral dipolar hard sphere model and we give 
some of its elementary properties. The subsection 2.1 is devoted to the definition 
of the model and to the computation of the multipole moments due to the two 
permanent dipoles moments ; all multipole moments are computed with the center 
of the sphere taken as origin. 

In subsection 2.2, we compute interaction energies between two enantiomers for 
few particular configurations ; in this subsection, a particular attention is pay 
to the chirodiastaltic interactions for colinear aligned total dipole moments and 
antiparallel dipole moment configurations. The discriminatory interaction ener gies 
are often computed by using the multipole expansion [1, [H, [13, H^, 41, IH, 47, 



l49l . 150( 1 ; thus, to facilitate the comparison with these works, the discriminatory 
interaction energies between chiral dipolar hard sphere have also been computed 
analytically by using the multipole expansion when the distance between both 
enantiomers is large in comparison to their diameter. However, if both enantiomers 
are at contact or if the distance between them is too small to use with accuracy the 
multipole expansion, then the interaction energies must be computed by summing 
all dipole-dipole interactions ; it is done so in subsection 2.2 for all configurations 
at contact and in section in all Monte-Carlo computations done with the Ewald 
method. 



2.1. Definition of the model and its multipole moments. 

The Chiral Dipolar Hard Sphere model consist of hard spheres of diameter a with 
two point dipoles Hi and /X2 located at points Oi and O2 inside the sphere. We 
define the molecular axis it of a chiral dipolar hard sphere as O1O2 = 2Lu with 
2L < a. The dipoles /^i and fj,2 are chosen such that ^Xa-u = and we set 

{/^i = ^ Ai 
fi2 = A/i (cos a fii + sin a {/jLi x ii)) 

where the parameters a G [— vr, vr] and < A < 1 define the second dipole fj,2 from 
the dipole /xi and the axis ii. We name Rectus (R-) bidipolar hard sphere (Rq-HS) 
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Figure 1. Representation the molecule-fixed local frame for the Sinister bidipolar hard sphere enantiomer 
(So-HS). The frame is defined from vectors u, fii and fX2 via Eas.l lU3l l ; | O1O2 |= 2L < a. The thick 
dashed vector is the projection of /X2 in the plane (Oi, e^, /Xj^). 

for < a < TT and Sinister (S-) bidipolar hard sphere (Sq,-HS) for — vr < a < 
(see figured]). In this nomenclature, it is clear that R(_q,)-HS = Sq,-HS ; molecules 
Rq,-HS and Sq-HS are enantiomers. 

Chiral parameters or chiral indices have been introduced to give a quantitative 



measure of the chirality of a molecule [20|, [SJ]- In particular, A.B. Harris and 



co-workers |54i | have defined chiral parameters as pseudoscalars constructed from 
the structure of molecules described via the representation theory of the three- 
dimensional rotation group 0(3) ; they have also related chiral parameters to the 
cholesteric pitch of generic liquid crystal models. For the chiral dipolar hard sphere 
model, we may defined a pseudoscalar from the three vectors ii, /xi and fi2 as 

X = Lu.{fj,2 X /^i) = LXfi^sina (2) 

and take x as the chiral parameter of this model ; we have xi^a — HS) = — x(Rq, — 
HS). If X = 0, then both enantiomers are equivalent and the bidipolar hard sphere 
model is achiral. This occurs if one of the following conditions is verified : (a) 
L = 0, then both point dipoles are located at the center of the sphere and the 
model reduces to the dipolar hard sphere model ; (b) A = 0, in this case, there 
is only one point dipole in the molecule and (c) sin a = then the dipoles and 
the molecular axis are all in the same plane which is a plane of symmetry for the 
molecule. 

A physical interpretation to | x* |= 1 is less obvious ; to obtain | x* |= 1 one needs 
to have A = 1, thus a symmetry is restored and then the chirality is perhaps less 
marked for these particular cases with A = 1 (for instance, see below subsection 
2.2 and Fig. 7). A trivial interpretation of values | x* |= 1 would be to say that for 
these value a maximum discrimination is achieved ; however, this is perhaps not 
fully correct. 

We define the molecule-fixed local frame as 

_ Ml + /^2 _ P 
\ Hl+ fj.2\ I -P I 



U 



(3) 
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Figure 2. Representation of the spherical components of the dipole and quadrupole moments as function of 
the chiral parameter x- The spherical components, represented as function of x, are multivalued because of 
the proprety x{'^) = xC"'"'*) (cf.Eq.l(2]l). Enantiomers R^-HS and Sa-HS have equal spherical components 
when an even symmetry occurs, that is to say for P and Q2i- One may note that for A = 1 one has Q21 = 0. 




-130 -so 



Figure 3. Representation of the spherical components of the dipole and quadrupole as function of a. As 
shown below, the knowledge of P, Q21 and Q22 as function of the parameters of the model determine fully 
all components of the spherical multipole tensor of an arbitrary order I (see below). 

where P is the total dipole of the molecule. In this molecule-fixed frame, the 
projections of the dipoles onto the frame axis are given by 



' /^i = fi{l + 2A cos a + A^) [—A sin a + (1 + A cos a) e^] 
/.i2 = + 2A cos Q + A^)""*^/^ [A sin a By + A(cos a + A) e^] (4) 
P = + 2Acosa + A^)!/^ 
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It is worthwhile to note that the total dipole of enantiomers Rq,-HS and Sq-HS 
in the molecule-fixed frame are the same, we have P(Ra-HS) = P(Sa-HS). One 
should note also that P = can be obtained only if A = 1 and a = vr, then x = 
and the model is achiral ; thus, all chiral bi-dipolar hard sphere have a non-zero 
net total dipole. 

The dipole distribution of the chiral hard sphere may be written as 



fi{r) = tii5{r + Lit) + tJ,2S{r - Lii) 



(5) 



and from multipole expansion for any dipolar distribution (see appendix A), the 
cartesian components of the quadrupole moment are given by 



2 V1 + 2A cosa + A2 



2Asina(A2-l) 
2AsinQ 
(A2 - 1) 



(6) 



and spherical components by 



Q22 = i 



Q21 
, Q20 




15 LXfisina 
27r Vl-F2Acosa-FA2 

15 Lfi{l - A^) 
Svr Vl + 2Acosa-hA2 



(7) 



For enantiomers, we have (522(Ra-HS) = — Q22fSra-HS) and (52i(R'a-HS) 
Q2i(Sa-HS). The magnitude of the quadrupole is [53] 



2_ -J r2, ,2 

2m I — fJ- 



4A2 sin2a + (l- A2)2n 
1 -h 2A COS a -h A2 



(8) 



Obviously, one has Q|(Ra-HS) = Q|(Sa-HS). In Eq.®, it might be surprising to 
find Q2 7^ for A = 0, since there is only a single dipole in the hard sphere in these 
cases. This stems from dependence on the choice of the origin in the computation 
of the quadrupole (the single dipole is shifted by —Lit from the center of the 
hard sphere - see also ref.[53], p. 69). This origin dependence highlight the fact 
that when setting A = 0, the chiral bidipolar hard sphere model do not reduce 
to the classical dipolar hard sphere model. On Figs.(l2][3]), we show the dipole and 
quadrupole components as functions of the chiral parameter x angle a ; on 
these representations, the influence of chirality is related to odd or even properties 
of the function. As shown on Figl2]l3l a difference between enantiomers Rq, and Sq, 
is found only in Q22- 
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The spherical components of the octopole are given by 



r Q33 = 



Q32 
Q31 = 



+ 2A cos a + A^) = \^L'Q^o 



1 105 

2 V Svr 



Q30 = -I J ^LV(1 + 2Acosa + X^)^/^ 
2 V 47r 



10 



(9) 



and its magnitude is 



33 



7 ^ 

m 



= ^ y: I Qsm \'= + 2Acosa + A^) = ^-^L^Qj, (10) 



It is apparent also on Eqs.([9] fT0|) that the octopole moments are the same for both 
enantiomers. More generally, by using Eq.(l5]) and the definition of the spherical 
component of multipoles Eq. ()A4p . we found 



il 



A sin a 



liLi^-^) Vl + 2Acosa + A2 



(l + (-l)')Y,i(e, 



(11) 



and, for < m < / — 1, we have 



Q 



Im 



/iL('-i) 



-^{l + m + l){l-m) 



(A^-(-l)0 + (l-(-l)')Acosa 
Vl + 2Acosa + A2 



-\-im- 



AsiuQ 



Vl + 2Acosa + A2 
The value of ^imif^x) may be easily computed, it comes 



(12) 



Y;r)i(^2 



21 + 1 



(-i)"^"'/- ";r'"';!;;""i:r for(i+,n) 

(/ — m)!!(t + m)!! 



even 



(13) 



for (/ + m) odd 



therefore, only one of the two contributions in the right handed side of Eq. (jl2p 
gives Qim- From this computation, we see that all spherical components can be 
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computed only from Qio, Q21 and Q22 ', more precisely, from Eqs. (|llllT2]) . we have 

2^ /Q22 
15 \fiL 



Qlm 




Ty = (l + (-!)') 



(i-(-iy) 



V^(/ + m + l)(/ -m) / Q21 
2 V /iL 



/{m+l)lea;J 



(14) 



/Q 



Since R(_fj)-HS = S^-HS, we have 



I2tt 



QlmiRa) - QUSa) = 2(1 + ("1) V^^'~'\/ TrQ22Y,„(e,.) 

15 



(15) 



thus, all spherical components of multipoles of enantiomers Rq-HS and Sq,-HS are 
equal, but those with / and m even and m 7^ 0. It is worthwhile to note that Qim 
is proportional to only one of the three components Qio, Q21 or Q22, because of 
Eq. p3|) . Eg . (1141) shows that the behavior of all components of multipoles are given 
(up to a numerical multiplicative factor) by Figs. ([2] -13]) as functions of the model 
parameters ; in particular the difference Qim(Ra) — QimiSa) is given in Figl2]||b), 
as function of x, and in FiglS^b), as function of a. 

As shown in Eq. (|14p all components of the multipoles dependent only on the three 
components Qio, Q21 or Q22, therefore it is advantageous to define the three re- 
duced geometrical parameters x*i P* and r* as 



X = - — K = A sm a 



P* = - = (1 + 2A cos a + A=^)i/2 



r* = =(1-A2) 
Lfj, 



(16) 



2.2. Interaction energies between two enantiomers. 

The interaction energy between two enantiomers can be computed in two ways. 
First, by summing up all interaction energies due to dipole-dipole interactions 
and/or, second, by using the multipole expansion. In the next section, we give some 
preliminary numerical results obtained with Monte Carlo Metropolis sampling of 
the phase space ; in these computations, all dipole-dipole interactions are taken 
into account explicitely and long ranged interactions are computed with the Ewald 
method [H,!!!]. 

In the present subsection, we give some analytical results for the interaction energy 
between two enantiomers by using the multipole expansion when r ^ a > 2L 
and we give also some analytical and numerical results obtained by computing 
explicitly all dipole-dipole interactions between two enantiomers at contact {i.e. 
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r = cr> 2L). 

From the multipole expansion, the interaction energy between the multipoles of 
orders li and I2 for two molecules of general shape is given by [13] 



»i> ^ ^ mi,m2,m 



ni,n2 



(17) 



where I = I1+I2 and Aij^ is a numerical constant, Qin the components of multipoles 
in the body-fixed frame, C{lil2l;rnim2m) the Clebsch-Gordan coefficients, -D^n 
are the rotation matrices, the Euler angles LOa = {(t>aOaXa) define the orientation of 
molecule a and uj = r. 

The first contribution to interaction energy is given by the total dipole-total dipole 
interaction and it is given by 



1 



f/ii(a = X„, 6 = Y„) = ^ [PaPb - 2,{Pa.r){Pb-r)] 



(18) 



where Xq, Yq, = Rq, or Sq/. This contribution is not discriminatory between enan- 
tiomers (see Eq.(j3])) and it is the longest ranged interaction between two chiral 
dipolar hard spheres. Therefore, many of the properties of this model will coin- 
cide with the properties of dipolar hard sphere systems that have been extensively 
studied by Jean- Jacques Weis in the last decade [5l|, |56|, |57|, |58j, |59|, |60j, ISlI, IS^, l63|, 

ei, m m, S S 0, El, Hi. 

The discriminatory interaction between enantiomers (Chirodiastaltic interaction 
H), may be defined by the difference A[/(^^) = f/(R„,R„) - [/(R„,S„) (where 

LU stands for Like — Unlike), we note also Aujj^^ = Uij^(Ra,^a) — Uhi20^cnSa) 
as the contribution of multipoles of orders li and I2 to the discriminatory chirodi- 
astaltic interaction. In the following, we choose the molecule at origin to be an Rq, 
enantiomer and all coordinates of vectors and positions are given in its molecule- 
fixed local frame given by Eq.(l3]) - this is equivalent to choose uji = (000) in 

Eq. (jl7p . From equation (llSp . we have Auj^j^^ 7^ only if I2 even, thus, the first 
non zero contribution to AU^^'^^ is given by the dipole-quadrupole interaction 
(/i = 1, I2 = 2) and after some algebra, it comes 



AU 



(LU) _ f Lfi 



12 I ^4 



with 



A sin a [15 M2(02te; ^0) + 12 Mi(02^2X2; ^0) + 9 Mo(02X2; ^)] 



(19) 



' M2((/>202X2; 0(1)) = sin^ e COS e [(1 + cos^ 62) sin 2^2 cos 2((/.2 
+2 cos 62 cos 2x2 sin 2{(f)2 — (/>)] 
Mi((/.2e2X2; 04>) = sin0(5cos2 ^ - 1) [cos 2^2 sin((/>2 - </>) 
-I- cos 62 sin 2x2 cos((/)2 - (p)] 

^ Mo (6^2X2; 6*) = (5cos3 6' - 3cos6')sin2 6'2sin2x2 



(20) 
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Figure 4. Representation of the two functions (21 and (22 as function of the reduced chiral parameter 
X* = x/^M^ s-Ki ^ functions of a. 



where L02 = {(P2G2X2) and lo = {(pO) are respectively the orientations of the second 
enantiomer and of the bond vector r in the molecule-fixed local frame of the first 
enantiomer. 

The quadrupole-quadrupole interaction contributes also to AC/^^^^ ; AC/22 
be written as 

Ao + 2OA1 + 30\/2^2 + 140yl3 + 35yl4 

(21) 

Bo + —Bi - 5V2B2 - —B3 



^^^^ --ll"^ 



f2i(A,a) 



with the functions f22 and f2i describing the dependence of AC/22 
parameters ; these two functions are given by 



(LU) 



on the model 



f22(A,a) 
f2i(A,a) 



sin^ a 



(l + 2Acosa + A2) 
(1 - A2)Asina 



p*2 

r*X* 



(l + 2Acosa + A2) P 



■*2 



(22) 



In Figsll](a,b), we represent these two functions versus the chiral parameter x and, 
in FigsH] (c,d) as functions of a. 
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(a) 



(6) 



Figure 5. Configurations of particular interest in dipolar system, (a) Colinear aligned total dipole mo- 
ments, these kind of configurations between two dipolar particle is responsible for chains formation at low 
temperature and low density : for this configuration, one has 9 = and 62 = 0. (b) Antiparallel total 
dipole moment : for this configuration, one has d = it/2 and 62 = tt. 

The functions A and B depend on orientations uj2 and lo. More precisely, we have 



' ^0 = (35 cos^ — 30 cos^ + 3) [(1 + cos^ 62) sin 2x2 sin 2(/)2 — 2 cos 62 cos 2x2 cos 2(/)2 
Ai = sin 6 (7 cos^ — 3 cos 6) sin 62 [cos 2x2 cos{(f)2 + (f) — cos 62 sin 2x2 sin((^2 + 0)] 
^2 = sin^ 6 (7 cos^ 6 — 1) sin 2(/) sin^ 02 sin 2x2 

^3 = sin^ cos sin 02 [cos 2x2 cos(<^2 — 3(^) — cos 02 sin 2x2 sm.{(j)2 — 30)] 
^ ^4 = sin^ 0[2 cos 02 cos 2x2 cos(2(/)2 — 40) — (1 + cos^ 02) sin 2x2 sin(202 - 



(23) 



and 



Bq = (35 cos^ — 30 cos^ + 3) sin 02 [cos 2x2 cos 02 + cos 02 cos 2x2 cos ^2] 
Bi = sin0 (7cos^ — 3 cos 0) [(1 + cos^ 02) sin 2x2 sin(202 - 0) 



+ 2 cos 02 COS 2X2 cos(202 — 0) 



2^2 



sin^ 02 sin 2x2 cos 1 



B2 = sir? (7 cos^ — 1) sin 02 [cos 2x2 sin(02 — 20) + cos 02 sin 2x2 sin(02 — 20)] 

i?3 = sin^ cos 0[(1 + cos^ 02) sin 2x2 sin(202 — 30) — 2 cos 02 cos 2x2 cos(202 — 30)] 

(24) 

Obviously, if the distance between two enantiomers is too small the multipole ex- 
pansion ceases to be accurate, however, the interaction energy between enantiomers 
may still be computed by summing all dipole-dipole interactions. The total inter- 
action energy between two enantiomers a and h has thus to be written as 



_g(a,6) ^ _g;(la,16) _^ _g;(la,26) _^ ^{2a,lh) _^ £;{2a,2fe) 



(25) 
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Figure 6. Interaction energies for two enantiomers in configurations given by Fig[5](a) at contact z = a. 
The angle (p is defined by cos = Ua-Ub- In all computations reported in these figures, /i = 1.0, L = 0.25o" 
and the curves correspond to an increase in the value of q by 15° ; the red lines correspond to a = 0°, 
45°, 90° and 135°, these values are indicated at edge of curves in figures. The energy < is the 



average energy over the orientation 11;, of enantiomers Y. (a) X and Y = Rq 
and A = 1.0 ; (c) X and r = Rq and A = 0.75 and (d) y = Sq and A = 0.75. 



and A = 1.0 ; (b) Y 



Table 1. Average energy over the orientation iib of 
enantiomers Y and location (pm of the orientation of iib- 
for which the minimun energy is reached, when Y = Rq 
for configurations given by Fig[5](a) at contact z — cr. 



a 


A = 1 


A = 0.75 


4>m 

A = 1 


A = 0.75 


0° 


-5.794 


-4.436 


0° 


0° 


15° 


-5.695 


-4.362 


11.2° 


10.4° 


30° 


-5.506 


-4.145 


22.4° 


20.8° 


45° 


-4.945 


-3.800 


33° 


31° 


60° 


-4.345 


-3.350 


43° 


40.2° 


75° 


-3.647 


-2.825 


51.6° 


48.4° 


90° 


-2.897 


-2.263 


59.2° 


55.4° 


105° 


-2.147 


-1.701 


65.8° 


61° 


120° 


-1.448 


-1.177 


71.4° 


65.2° 


135° 


-0.848 


-0.727 


76.6° 


67° 


150° 


-0.388 


-0.381 


81.2° 


62.8° 


165° 


-0.098 


-0.164 


85.6° 


11° 



where denotes the interaction energy between the dipole i of enantiomer a 

with dipole j of enantiomer 6 ; it is given by 



j^(ia,jb) 



1 



r.«J 13 



ab 



(26) 



and the distance vectors between point dipoles are 
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Even for simple relative configurations, as the ones shown on Figl5]^a-b), the total 
interaction energy between two enantiomers is a bit complicated when computed 
by summing all dipole-dipole interactions. For instance, the interaction energy 
between two Rq, enantiomers in a configuration given by Figl5][^a) is 



^{a,b) 



(l + 2(f )(l-COS,/.))5/2 



p*'2 _j_ 



P 



*2 
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P 



*2 



(l + COS(^)-l 



+6 1 -]x*sm<P-2(^ I I cos( 



L 



(2 + 2 COS (/) + sin^ 



(28) 

where cos(j) = Ua-Uh and (x*,-P*,r*) given by Eq. (ll6p . The chirodiastaltic interac- 
tion between enantiomers Rq, and Sq,, computed by summing the four dipole-dipole 
interaction energies, for the configuration given on Figllja) is 



^E^LU) ^ ^(R„,R„) _ ^(R„,S„) 



(l + 2(f) (l-COS,/.))5/2 
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Figure 7. Interaction energies for two enantiomers in configurations given by Fig[5](b) at contact r = a for 
A = 1. In all computation, = 1.0, L = 0.25(T and notations arc the same as in Fig[6] For the configuration 
given by Fig[5](b), the angle is defined by cosc^ = iia-r and the angle 02 by cos 02 = Ua-Ub- (^-b) : 
= 0° and (c-d) : = 45°. 




Figure 8. Interaction energies between two enantiomers in configurations given by Fig[5] (b) at contact 
r = cr for A = 0.75. Notations are the same as for Fig[7] 
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For \ z \^ L, a. Taylor expansion of Eq. (j29p gives 



z 



2 






(30) 



AU^;''>{92 = 0,e = 0) + AU^^''\62 = o,e = o) + 



o 




) 



where we have set (/> = i;^2 + X2 in Eqs. (23-24) to obtain AC/22 • Substracting 
Eq. (j30p from Eq. (129p . may allow to check the degree of accm'acy of the multipole 
expansion for the configuration of Figl5] (a). 

On Figsini we represent interaction energies for enantiomers in the colinear aligned 
total dipole moment configuration (FigHJ-a). These energies have been computed by 
summing all dipole interactions between the two enantiomers. To allow simple com- 
parison between interaction energies, for different values of a, we have substracted 
to all curves the average, over the orientation iih, of the interaction energy. For 
all A, a minimum in the interaction energy between two different enantiomers is 
always obtained for (j) = (iia and Uh parallel) ; the anti-parallel configuration 
of molecular axis iia and iify {(j) = vr) is a relative minimum for the interaction 
energy for all A < 1 (see Figj6)^b-d) - if A = 1, the configurations for (p = and 
(j) = TT have the same energy) ; one should note also that these minima are relative, 
because the configurations like in Figl5] (a) do not permit to reach the smallest 
distance between two point dipoles. The minimum interaction energy between two 
enantiomers of same kind is not obtained for parallel or anti-parallel configurations 
of the molecular axis. The value 4>rn of the angle for which a minimum interaction 
energy is obtained, depends on parameters A and a (see Figl6]^a-c)), several values 
are given in Table [1] ; average values of < E-^^ >^ are also reported in this Table. 
Therefore, in homochiral systems (systems that are composed with exclusively one 
enantiomer) with a large enough total dipole moment and in thermodynamical 
conditions where dipolar systems form chains (sil. [Hs!]. we may observe formation 
of chains too, and the molecular axis of the enantiomers may have helical arrange- 
ment. Addition of chiral enantiomers in the system will reduce the range of the 
helical arrangement in the chains. For A = 1, the minimum at (j) = (j)m has the 
same value than the one found at = (pm — tt, therefore right and left-handed 
helical arrangement of the molecular axis in the chains will compete in homochiral 
systems. 

For the configuration given in Figl5^b) and for r ^ L (9 = Tr/2 and 62 = it), the 
multipole-multipole chirodiastaltic interaction is found as 



On Figsl7]and[8l we represent interaction energies for enantiomers in the antipar- 
allel total dipole configurations (Figl5]-b) for = and it /A. On Figsd (a,c) and 
[8](a,c), we represent interactions energies between two enantiomers for A = 1 
and 0.75 and on Figsl7](b,d) and[8](b,d) we represent interactions energies between 






(31) 
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enantiomers Rq and Sq,. The coupling between tj) and (j)2 renders quite complicated 
the energy landscape of antiparallel configurations. The relative minima for an- 
tiparallel configurations (Figl5]-b) are smaller than the minima found for colinear 
configurations (Figl5]-a), therefore the colinear and antiparallel configurations are 
in drastic competition by comparison to dipolar hard sphere systems. 
For two identical enantiomers, when (j) = n7r/2 (with n = —1,0, 1 or 2), the loca- 
tion of minimum interaction energies are independent on a. For instance, if = 0, 
the relative minima for two Rq enantiomers are always obtained for 02 = or vr ; 
if A = 1, these two minima have the same value, but if A < 1 the two minima are 
not equivalent and their value dependent on a. On FigiH we give the value of the 
minimum energy for </) = and (/)2 = or vr for A = 0.75 and 0.5 ; as shown on this 
figure for values of a > am(A) the configuration (/>2 = vr is more stable than the 
configuration (/>2 = 0. The amplitude of variation of interaction energies with (j)2 
are larger for antiparallel configurations than variation with (p for aligned colinear 
configurations. 

The locations of relative minimum interaction energies between enantiomers Rq 
and Sq, depend not only on (f) and (j)2-, but also on a. For instance on FiglT^b) 
{(j) = 0) and for a = 57r/12 (75°), we found four different minima located at 
(t)2 ^ ±141° and ±39°. 

The computations of interaction energies between two chiral dipolar hard spheres 
show that a lot of different configurations may be in competition, even for sim- 
ple configurations as the ones shown on FigsiSl In the next section, we give some 
preliminary results obtained with Monte Carlo simulations for several cases. 



3. Preliminary Monte Carlo results 

As shown in the previous section, interaction energies between two enantiomers as 
function of their relative orientation may have a lot of relative minima ; therefore, 
despite the apparent simplicity of the model, the thermodynamical properties chi- 
ral dipolar hard sphere systems are quite complicated. In particular, because the 
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energy landscape of the interaction between two enantiomers shows large varia- 
tions, one may predict that frustration phenomena will play an important role in 
phases at low densities, low temperatures and also in solid phases. 
We define the reduced quantities as follow : for volume of the simulation box 
V* = V/a^, density p* = pa^ = Na^/V, pressure p* = pa'^/kT and dipole mo- 
ment p* = (p/kTa^)^^"^ . For a given model of enantiomer, defined by (A, a) and L, 
the reduced total dipole moment P* of the molecule is related to p* according to 
Eq.([4]). In all computations, we have chosen L* = L/a = 1/4 as in the numerical 
computations done in the previous section (see FigslMH]). At the end of this section, 
a short discussion on the infiuence of L on structures is presented ; this influence 
is illustrated on two cases : L* = 0.025 and L* = 0.425. For notational convenience 
the asterisks will be drop in the following. 

Mixtures of enantiomers Ra-HS and So-HS are also of interest, we define densities 
/9r = N^/V and ps = N^/V respectively for enantiomers Rq-HS and Sq,-HS and 
the ratio X = Ps/Pr = N^/N^. The ratio X is related to the enantiomeric excess 
(ee) in enantiomer Rq-HS by ee = 1 — X. In the following, we have considered only 
two mixtures of enantiomers : homochiral R^ systems {X = 0) and racemic Ro-Sq 
systems {X = 1). 

Monte Carlo simulations have been performed in canonical (NVT) and isobaric 
(NPT) ensembles with system size N = + Ng = 512 or 1000. Periodic bound- 
ary conditions are used and dipolar interaction energies between two enantiomers 
is computed with Eqs. (j25ll271) where long ranged contributions are taken into ac- 



count by using the Ewald method [5lf|. 



From a technical point of view, as a consequence of the numerous number of rela- 
tive minima in interaction energies, the convergence of the Monte Carlo algorithm 
and the sampling of the phase space will require a large number of MC cycles (one 
MC cycle corresponds, on average, to one trial move per particle for NVT simu- 
lations ; for NPT simulation, an additional trial move of the volme of the box is 
done every MC-cycle). For most results presented in this section, between 5 x 10^ 
and 10^ MC cycles have been achieved for the relaxation and convergence from an 
initial condition and averages have been accumulated over about 10^ and 2 x 10^ 
MC cycles, respectively for systems with 1000 and 512 particles. 
Periodic boundary conditions are used and interaction energies are computed with 
the Ewald method [HI] by summing all dipole-dipole interactions between enan- 
tiomers. The average energy U is computed as 



C/=(y^("'''M (32) 




where E^'^'^^ is the interaction energy between enantiomers a and h Eqs. (j25ll27|) and 
< . > is the average of the MC sampling. Some computation have been done with 
an external electric field Eq = EqBz, the energy in the external field is computed 
as 



V=(^-^P^'^\Eo'j = l^-^U^^ + f^i^^Eo"^ (33) 
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Tabic 2. Average energies and eigenvalues S4. , So and S_ for homochiral-Ro, and racemic systems in external electric field E — Ee~ 
{E — 10) obtained in canonical (NVT) Monte Carlo simulations. The numbers in brackets give the accuracy on the last digit of the 
averages. For all systems A^^ — 512, /j. — 1.0 and p — 0.05 (or p — 0.5). The parameters a and A define the molecular structure of 
enantiomers as given by Eg. (fit . 
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For these computations, the density is p = 0.5. 




Figure 10. Representation of gooo('') for several value of the pressure obtained by Monte-Carlo sampling 
of the isobaric (NPT) ensemble of homochiral-Ra systems for N = 1000, fj, = 0.5, a = 60° and A = 1.0. 



Possible orientational order can be established by computing the eigenvalues of the 
tensor 
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Figure 11. Snapshots of homochiral-Rc and racemic systems in external electric field. For all systems 
= 512, /i = 1.0, p = 0.05, A = 1.0 and E = Ee^ {E = 10) ; the solid vertical line indicates the 
direction of the field, (a-b); homochiral-Ra with a = 35° ; (c-d): racemic Rc-Sa system with a = 35° and 
(e-f)homochiral-RQ with a = 150°. Ra-HS enantiomers are represented in blue and Sa-HS enantiomers 
in red. In (b), (d) and (f), we represent, for the same configurations, molecular axis and segment in the 
direction of dipoles for each enantiomer ; dipoles of Rq-HS enantiomers are represented in blue and dipoles 
of Sq-HS enantiomers in red. 

where is the cartesian component of a unit vector x of molecule i. The largest 
eigenvalue S+ of Qa/s is used as an order parameter for the isotropic-nematic phase 
transition in liquid crystals. It is also worthwhile to note that if, all unit vectors x 
are perpendicular to a given fixed direction then, the eigenvalues are 5+ = = 1/4 
and = —1/2. In the present work, we have restricted the computation of Qajs 



April 4, 2009 2:11 Molecular Physics Mazars'chiral'condmat 



20 

Tabic 3. Preliminary results for Montc-Carlo computations in 
the isobaric [N PT) ensemble of homochiral-Ro^ systems. These 
computations are performed with N — 1000, f-L — 0.5, a — 60° 
and A — 1.0. Notations are the same as in Tabic [2] p is the 
reduced pressure and < ^9 > is the average value of the density. 
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Figure 12. Snapshots of chiral dipolar hard sphere systems for A'^ = 1000, fi = 2.0, p = 0.1, a = 120° 
and A = 1.0. The notations are the same as those in Figs llll (a-b); homochiral-Ra system ; (c-d): raccmic 
Ra-Sa system. In (a) and (c) hard spheres are represented as balls and molecular axis O1O2 by a black 
segment. No external electric field are applied to the systems. 



for the molecular axis u of enantiomers. 

To study the structure of the systems we have computed the angular pro jections 
on rotational invariants of the molecular pair distribution function (sol. l6ll. [tsI. 74] 
; in the present work, we have computed only gooo{r), hiiQ{r) and hii2{r) for two 
sets of points and vectors. The first set is the molecule center and the orientation of 
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1.0 2.0 3.0 4.0 5.0 0.5 1.0 1.5 2.0 2.5 3.0 



r/a r/o 

Figure 13. Representation of projections of molecular pair distribution functions computed by Monte 
Carlo sampling of the {NVT) ensemble of racemic Rq-Sq systems. In these computation, A'^ = 1000, 
fi = 1.5, p = 0.1 and a = 30° ; (a-b) A = 1.0 and (c-d) A = 0.75. The functions represented in (a) and (c) 
are computed for the set (O, P) and the functions represented in (b) and (d) for the set (Oi,/li). In all 
figures, pair distributions functions for Like-Like pairs of molecules are represented in black and those for 
Like- Unlike pairs in red. 

the total dipole moment : (O, P) ; the second set is the site Oi and the orientation 
of the dipole /xi : (Oi, p,i). To study the influence of the chirodiastaltic interaction 
on structure of mixtures of both enantiomers, it will be necessary to compute pro- 
jections of the pair distribution functions on rotational invariants of higher order 
since the discriminatory interaction involves dipole-quadrupole interactions. 
Table [21 summarizes average energies and eigenvalues of Qa/3 Eq. fjM]) for 
homochiral-Rc and racemic systems in external electric field E = Ez (E = 10), 
for several values of a and A. On Table [3l we report some preliminary results of 
Monte-Carlo computations in the isobaric (NPT) ensemble of homochiral-Ra sys- 
tems and on FigllOl we show gooo{r) obtained with these computations. 
On FigdU we give few snapshots of homochiral-Ro and racemic systems in an 
external field. The results reported on table [2] show that energies and the nematic 
order parameters for the molecular axis are almost independent on the composition 
of systems. The eigenvalues 5+, So and show that on average the molecular axis 
of enantiomers are perpendicular to the direction of the electric field, this results 
is quite obvious since the total dipole moment, as defined in Eqs. ()3|4p . is per- 
pendicular to the molecular axis this tendency decreases as a increases, since the 
total dipole moment decreases. As shown on the snapshots in Fig. ijlip . li a < 60° 
molecules tend to condense in a large columnar cluster, for both homochiral-Ra 
and racemic systems in an external electric field (cf. Figllll (a-d)) ; if a > 70° the 
columnar cluster does not form (cf. Figllll (e-f)). There are at least two facts that 
may help to understand qualitatively the disappearance of the columnar clusters 
as a increases, while A and /i are fixed. First, as a increases, the total dipole of 
each molecule decreases and thus formations of chains in the direction of the ex- 
ternal field is less favoured. Second, for large enough a, configurations with two 
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Figure 14. Snapshots of chiral dipolar hard sphere systems for homochiral-Ra system with A'^ = 512, 
= 1.5, p = 0.05, a = 30° and A = 1.0 and for two value of L. The notations are the same as those in 
Figs llll (a-b): 2L = O.OSct ; (c-d): 2L = 0.95(t. In (a) and (c) hard spheres are represented as balls and 
molecular axis O1O2 by a black segment. No external electric field are applied to the systems. 



dipoles belonging to two different molecules may be more easily found, on average, 
in an antiparallel configuration in a plane perpendicular to the field, while the total 
dipole of both molecules can still be parallel to the field. At present, it is not clear 
how the columnar clusters form and which mecanism favour their formation ; fur- 
ther investigations are needed to be conclusive. Homochiral-Ro and racemic Rq-Sq, 
systems in an external electric field exhibit the same properties, however, a close 
inspection of the snapshots of FiglTT] (a-d) shows that the internal structure of the 
columnar clusters is not the same in homochiral and racemic systems. In homochi- 
ral systems, the molecules in the columnar cluster have an helical arrangement 
(in agreement with the results shown on FigjG]) ; in racemic systems, these helical 
arrangement in the columnar clusters are less marked because of the Like-Unlike 
interactions. To obtain more quantitative descriptions of the internal structure of 
the columnar clusters, some particular order parameters and distribution functions 
have to be defined. 

If there is no external electric field, columnar clusters are not formed instead some 
complicated structures are formed at low density and large enough // ; this has oc- 
cured for all values of a that have been considered in these prelyminary results. On 
Fig ll2l snapshots of homochiral-Ro, and racemic Rq-Sq, systems are represented, for 
these systems a = 120° and no external electric fileds are applied. Again, in these 
computations homochiral and racemic systems behave very similarly, however, an 
inspection of the pair distribution functions show that there is small differences 
between pair distributions functions for Like-Like pairs of molecules and those for 
Like-Unlike pairs (see for instance Fig llSp . The comparison between pair distribu- 
tion functions computed for A = 1.0 and 0.75 shows that this parameter have an 
important influence on the local structure (cf. Fig llSp . In particular, when A ^ 1.0, 
interactions between the dipoles fii of two enantiomers are favoured and so an 
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asymmetric arrangement of the molecular axis. 

There are three important parameters in this chiral model : A, a and L. So far in 
this section in the simulations reported (and also in numerical computations done 
in subsection 2.2), we have set always Lja = 1/4 and showed how energy minima 
and structures depend on A and a. They will depend also on the ratio L/a : for 
instance, if L/cr is small, chains with colinear aligned total dipole moment will be 
favoured, while, if L/a is close to 1, structures with colinear aligned molecular axis 
will be favoured, this behaviour is illustrated on Figll41 



4. Perspectives 

A full determination of the structures and clusters at low densities and temper- 
atures will require a larger sampling of the phase space ; this can be achieved 



with longer runs and by implementing trial cluster moves [75|, |76|, |77|]. Nematic 
and helical-cholesteric order parameter and also the bond orientational order 
parameters [t^ will be useful to investigate the properties of structure at low 
density as function of the enantiomeric excess and model parameters, homochiral 
systems have interest by themselves. Same order parameters and simulations in 
the isobaric ensemble (NPT) will serve to study the crystal phases of homochiral 
systems and also as functions of ee. Works and computations in these directions 
are already in progress. 

Another chiral model may be defined from the computation of the multipole 
moments done in section 2, and in agreement with earlier studies on chiral 
discrimination due to multipole-multipole interactions i, H Ei]. Accord ing to 
Eqs. (4,7), the three components, Qio, Q21 and Q22 are independent, thus an hard 
sphere model with dipole-dipole, dipole-quadrupole and quadrupole-quadrupole 
interactions [tI, HO] can be built by using these three multipole components as 
parameter for the model instead of (A,a,L). Such model will be slightly different 
from the bi-dipolar hard sphere model (and also from the four charges model), 
since higher order multipole-multipole interactions will not be taken into account 
in this model. Such systems has not been studied in the present work. 
With the bidipolar model, one may easily build a model of chiral liquid crystal by 
changing the short ranged steric repulsion by replacing the hard sphere potential 

or ellipsoid hard potential (s^ . [83 |. or even by a 



by a spherocylinder [8ll . 
soft Gay-Berne potential j85l. l86l|. For these chiral liquid crystal models with two 
dipoles, all the results of section 2 hold and in particular : the computations of 
the multipole components and interaction energies between to enantiomers. The 
modifications of the phase diagram of these systems with the introduction of these 
chirodiastaltic interactions can be of particular interest for the physics of liquid 
crystals. 
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Appendix A. Multipole expansion for a dipole distribution. 

In this appendix, we give general formulas for the multipole expansion of a dipole 
distribution /x(t") located in a finite region 0, of the space. From classical elec- 



tromagnetism theory [521 . l53l | and choosing the origin inside the electrostatic 



potential $(r), due to the dipole distribution in Q, is given by 



(Al) 



Since we are interested in the potential at r outside the dipole distribution, we set 
r > r' and r' and f, the unitary vectors that define directions of r' and r. The 
relation 

, _ I - E E (^) ;:iffiTY,.»(i-')YF„(i-) (A2) 

I 1=0 m=-l ^ ^ 



is reported in Eq. flAlh and the multipole expansion for a dipole distribution is 
given by 



- E E (^) (A3) 

;_n ^ — _i \ / 



1=0 m=-l 



with the m*'^ component of the spherical multipole moment tensor of order / 
defined by 



Qim = [ f^{r').Vr' \r" Yimir')] dr' (A4) 



Eqs. (|A3tiA4p are the multipole expansion for a dipole distribution. In the following, 
we give explicitly the Qim up to Z = 3 for any dipole distribution. From Eq. (|A4p . 
and since any dipole distribution can also be represented as a real charge density, 
it is easy to show that 

Ql,.m = i-irQlm (A5) 



A.l. The total charge in fl : Qqq. 
From Eq. ()A4p . we have 

Qoo= / /i(r).V,[Yoo(f)]cir = (A6) 
Jn 



Of course, this result is obvious for any dipole distribution. 
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A. 2. The dipole moment tensor : Qim- 

The components Qn and Qiq are given by 




r sm f/e 



i4> 



dr 



Qio = \j^(yj M(r-).Vr [rcose]dr 



thus, we find 




Qii = -\^[Px + iPy 



Q 



10 



47r 



27 



(A7) 



(A8) 



Where P is the total dipole in Q given by 

P = ^i{r)dr 
Jn 

This result is exactly the same as the one obtained for a charge distribution p{r) 
with the total dipole defined by 

P = rp{r)dr 
Jn 



A. 3. The quadrupole moment : Q^m- 

The same derivation gives for the spherical components of the quadrupole 



1 /15 
2V 2^ 



Q22 

g2i 



Q20 — \ -r 





/ xpxdr - I ynydr + i / {ypx + xpy)dr 
Jn Jn Jn 

/ (x/i^ + zij.r,)dr + i / {ypz + z^iy)dr 
Jn Jn 

2 / zfXzdr — / xfixdr — / y^ydr 
Jn Jn Jn 



(A9) 



The relations between spherical and cartesian components are given by 53l |. 



Q22 = \l {qxx - qyy + '^iqxy) 



Q 



21 




{Qxz + iqyz) 



(AlO) 
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and the traceless quadrupole moment tensor in cartesian coordinates can be written 
as 



(All) 



A. 4. The octopole moment : Qsm- 

The spherical components of the octopole moment may also be computed from 
Eq. (|A4p ; after some algebra, we found 



Q 



33 



3 /35 

4 V 47r 



Q 



1 105 



32 



Q 



2 V Svr 



1 21 



31 




8 V 47r 



/ {fi^r: + iiJ.y){x + iyfdr 
Jn 

/ ^jLz{x + iyf'dr + 2 / {fi^ + iiJ,y)z{x + iy)dr 
Jn Jn 

8 / ^izz{x + iy)dr + / {^^ - ifj,y){x - iyfdr (A12) 
Jn Jn 



+ 4 / {fi^iz^ - x^) + i^iy{z' - y^))dr 
Jn 



/n 2 2 2 

Hzi^z^ -x-y'' 



)dr -2 I zifx^x + fiyy)dr 
Jn 



and the cartesian components of the octopole moment are given by 

If 5 
^ 3 

-2{r.tJ,){raSf3^ + r/3(5a^ + r^S^a)] 



(A13) 



